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Abstract 

We complement the theory of tick-by-tick dynamics of financial markets based on a 
continuous-time random walk (CTRW) model recently proposed by Scalas et al [4], 
and we point out its consistency with the behaviour observed in the waiting-time 
distribution for BUND future prices traded at LIFFE, London. 
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1 Introduction 

In financial markets, not only prices can be modelled as random variables, but 
also waiting times between two consecutive transactions vary in a stochastic 
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fashion. This fact is well known in financial research. In his 1973 paper [1], 
Peter Clark wrote: "Instead of indexing [time] by the integers 0,1,2,. . ., the 
[price] process could be indexed by a set of numbers ti, t 2 , t 3 , . . ., where these 
numbers are themselves a realization of a stochastic process (with positive 
increments, so that t\ < t 2 < t% < . . .)." 

Till today, there have been various studies on the nature of the stochastic 
process generating the sequence of the tj. In Clark's approach, the variable 
t is not a physical time but an economic variable, the trading volume, an 
observable whose increments represent the market intensity. 

Lefol and Mercier have written a review [2] on the works inspired by Clark's 
seminal paper. In a review by Cont [3], readers can find pointers to the relevant 
literature and a description of the main research trends in this field. 

In a recent paper, Scalas et al. [4] have argued that the continuous time random 
walk (CTRW) model, formerly introduced in Statistical Mechanics by Montroll 
and Weiss [5] (on which the reader can find further information in Refs. [6- 
12]), can provide a phenomenological description of tick-by-tick dynamics in 
financial markets. Here, we give further theoretical arguments and test the 
theoretical predictions on the waiting-time distribution against empirical data. 

The paper is divided as follows. Section 2 is devoted to the discussion of a 
new form for the general master equation in the case of non-local and non- 
Markovian processes. In Section 3, the conditions for the derivation of the 
time-fractional master equation are given. The Mittag-Lefner function plays 
a central role in this respect. In Section 4, the theoretical predictions on the 
waiting-time distribution are compared to market data: high-frequency BUND 
future prices traded at LIFFE0 in 1997. Finally, the main conclusions are 
drawn in Section 5. 



2 The general master equation and the "memory function 



Throughout this paper the variable x represents the log-price. In other words, 
if S is the price of an asset, x = log 5*. The reason for this choice is explained 
by Scalas et al. [4]; it is essentially due to the fact that, rather than prices, 
returns are the relevant variable in finance. The physicist will recognize in x 
the position of a random walker jumping in one dimension. In the following, 
we shall often use the random walk language. 

Let us consider the time series {x(ti)} , i = 1,2, ... , which is characterised 
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by </?(£, r), the joint probability density of jumps & = x(tj) — x(tj_i) and 
of waiting times = tj — The joint density satisfies the normalization 
condition J °° <p{£,t) d£ dr — 1 . Relevant quantities are the two prob- 
ability density functions (prf/'s) defined as A(£) := ip(^,r) dr , ip{r) : = 
f^™ T ) d£ , and called jump pd/ and waiting-time pdf, respectively. 

The CTRW is generally defined through the requirement that the 7* are 
identically distributed independent (i.i.d.) random variables. Furthermore, in 
the following we shall assume that the jump pdf A(£) is independent of the 
waiting-time pdf, ipij) , so that the jumps & (at instants tj , i = 1,2,3,...) are 
i.i.d. random variables, all having the same probability density A(£) . Then, 
we have the factorization </?(£, r) = A(£) ipij) . For convenience we set t = . 

The jump pd/ A(£) represents the pdf for transition of the walker from a point 
i to a point x + £ , so it is also called the transition pdf . The waiting-time pdf 
represents the pdf that a step is taken at the instant + t after the previous 
one that happened at the instant ti-i , so it is also called the pausing-time pdf . 
Therefore, the probability that r < tj — < r + dr is equal to ^(r) rfr . 

The probability that a given interstep interval is greater or equal to r will be 
denoted by ^(r) , which is defined in terms of ipij) by 

OO T 

*(r) = / m dt' = l-J 4>(t') dt' , V(r) = • (2-1) 

r 

We note that Jq ip(t') dt' represents the probability that at least one step is 
taken at some instant in the interval [0, r), hence ^(r) is the probability that 
the diffusing quantity x does not change value during the time interval of 
duration r after a jump. We also note, recalling that t — , that ty(t) is the 
survival probability until time instant t at the initial position xo = . 

Let us now denote by p(x, t) the pdf of finding the random walker at the 
position x at time instant t . As usual we assume the initial condition p(x, 0) = 
5(x) , meaning that the walker is initially at the origin x — . We look for the 
evolution equation for p(x, t) , that we shall call the master equation of the 
CTRW. Montroll and Weiss [5] have shown that the Fourier-Laplace transform 
of p(x, t) satisfies a characteristic equation, now called the Montroll- Weiss 
equation. However, only based upon the previous probabilistic arguments and 
without detour onto the Fourier-Laplace domain, we can write, directly in the 
space-time domain, the required master equation, which reads 

t r +oo 

p(x,t) = 6(x)V(t) + J ip(t-t') J \(x- x')p(x',t')dx' dt'. (2.2) 

L-oo 

The spatially discrete analogue of this purely integral form of the master 
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equation is quoted in Klafter et al. [13] (see also Ref. [14]). We recognize from 
Eq. (2.2) the role of the survival probability ^(t) and of the pdf's ip(t) , X(x) . 
The first term in the RHS of (2.2) expresses the persistence (whose strength 
decreases with increasing time) of the initial position x = 0. The second 
term (a spatio-temporal convolution) gives the contribution to p(x, t) from 
the walker sitting in point x' G R at instant t' < t jumping to point x just at 
instant t , after stopping (or waiting) time t — t' . 

Now, passing to the Fourier-Laplace domain, we can promptly derive the cele- 
brated Montroll- Weiss equation [5]. In fact, by adopting the following standard 
notation for the generic Fourier and Laplace transforms: 

+oo oo 

F {/(*);«} = /(«)= / e tKX f(x)dx, L{/(0; s} = /(«) = / e~ st f(t) dt , 

-oo 

we get from (2.2) the Montroll- Weiss equation: 



Hereafter we present an alternative form to Eq. (2.2) which involves the first 
time derivative of p(x, t) (along with an additional auxiliary function) so that 
the resulting equation can be interpreted as an evolution equation of Fokker- 
Planck-Kolmogorov type. 



For our purposes we re-write Eq. (2.3) as 



Ms) 



sp(K, s) 



A(«) 



P(K,S) 



where 



l-tp(s) *(s) 



s i/j(s) 

Then our master equation reads 



ip(s) 1 - sV(s) 



(2.4) 
(2.5) 



J $(t-t')-^p(x,t')dt' = -p(x,t)+ J X(x-x')p(x',t)dx' , (2.6) 



where the "auxiliary" function $(t) , being defined through its Laplace trans- 
form in Eq. (2.5), is such that ^(t) = / * $(t - t') ip(t') dt' . We remind the 
reader that Eq. (2.6), combined with the initial condition p(x, 0) = S(x) , is 
equivalent to Eq. (2.4), and then its solution represents the Green function or 
the fundamental solution of the Cauchy problem. 

From Eq. (2.6) we recognize the role of "memory function" for <&(£) . As a 
consequence, the CTRW turns out to be in general a non-Markovian process. 
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However, the process is " memoryless" , namely "Markovian" if (and only if) 
the above memory function degenerates into a delta function (multiplied by 
a certain positive constant) so that \P(t) and ip{t) may differ only by a multi- 
plying positive constant. By appropriate choice of the unit of time we assume 
$(s) = 1 , so $(t) = S(t) , t > . In this case we derive 

$(s) = *(s) = -±-, so V(t) = tf(t)=e-*, t>0. (2.7) 
Then Eq. (2.6) reduces to 

— p(x, t) = -p(x, t) + J \(x- x')p(x',t)dx' , p(x,0)=5(x). (2.8) 

— oo 

This is, up to a change of the unit of time (which means multiplication of 
the R.H.S by a positive constant), the most general master equation for a 
Markovian CTRW; it is called the Kolmogorov-Feller equation in Ref. [15]. 

We note that the form (2.6), by exhibiting a weighted first-time derivative, 
is original as far as we know; it allows us to characterize in a natural way a 
peculiar class of non-Markovian processes, as shown in the next Section. 



3 The time-fractional master equation for "long-memory" pro- 
cesses 



Let us now consider " long- memory" processes, namely non-Markovian pro- 
cesses characterized by a memory function <3>(t) exhibiting a power-law time 
decay. To this purpose a natural choice is 

^ = T{ llpy 0</3<l. (3.1) 

Thus, $(£) is a weakly singular function that, in the limiting case /3 — 1 , 
reduces to <3>(i) = S(t) , according to the formal representation of the Dirac 
generalized function, S(t) = r7r(0) , t > (see e.g. REf. [16]). 

As a consequence of the choice Eq. (3.1), we recognize that (in this peculiar 
non-Markovian situation) our master equation (2.6) contains a time fractional 
derivative. In fact, by inserting into Eq. (2.4) the Laplace transform of $(£) , 
$(s) = l/s 1_/3 , we get 

s^p(k,s) - s^ 1 = \\(k) - ll p(k,s) , < /? < 1 , (3.2) 
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so that the resulting Eq. (2.6) can be written as 

— p(x, t) = -p(x, t) + J \{x - x')p(x',t)alx' , p(x,0) = 5(x), (3.3) 

— oo 

where d^/dt^ is the pseudo-differential operator explicitly defined in the Ap- 
pendix, that we call the Caputo fractional derivative of order (3 . Thus Eq. 
(3.3) can be considered as the time-fractional generalization of Eq. (2.8) and 
consequently can be called the time-fractional Kolmogorov- Feller equation. 
We note that this derivation differs from the one presented in Ref. [4] and ref- 
erences therein, in that here we have pointed out the role of the long-memory 
processes rather than that of scaling behaviour in the hydrodynamic limit. 
Furthermore here the Caputo fractional derivative appears in a natural way 
without use of the Riemann-Liouville fractional derivative. 

Our choice for implies peculiar forms for the functions ^(t) and ip(t) 
that generalize the exponential behaviour (2.7) of the Markovian case. In fact, 
working in the Laplace domain we get from (2.5) and (3.1) 

from which by inversion we obtain for t > 

= Ep(-tf>) , m = -j t M~t P ), 0<(3<l, (3.5) 

where Ep denotes an entire transcendental function, known as the Mittag- 
Leffler function of order /3 , defined in the complex plane by the power series 

oo n 

For detailed information on the Mittag-Leffler-type functions and their Laplace 
transforms the reader may consult e.g. [17-20]. We note that for < (3 < 1 
and 1 < f3 < 2 the function ^(t) appears in certain relaxation and oscillation 
processes, then called fractional relaxation and fractional oscillation processes, 
respectively (see e.g. Refs. [18,19,21,22] and references therein). 

Hereafter, we find it convenient to summarize the features of the functions 
ty(t) and ip{t) most relevant for our purposes. We begin to quote their series 
expansions and asymptotic representations: 

oo jjin 

sin(/3vr) T{f3) 



7T t? 



t — > oo , 
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and 

^S (_1)n r(/5n + /5)' *-° m 



sin(/?7r) + 1) 

t — > oo . 



71 t^ 1 ' 

The expression for tp(t) can be shown to be equivalent to that one obtained in 
Ref. [14] in terms of the generalized Mittag-Leffler function in two parameters. 

In the limit for f3 — > 1 we recover the exponential functions of the Markovian 
case. We note that for < (3 < 1 both functions ip(t), ^(t), even if losing 
their exponential decay by exhibiting power-law tails for large times, keep the 
"completely monotonic" character. Complete monotonicity of the functions 
ip(t), V(t), t > 0, means: 

(-l) n ^^)>0, (-l) n ^m>0, n = 0,l,2,... (3.9) 

or equivalently, their representability as (real) Laplace transforms of non- 
negative functions. In fact it can be shown for < (3 < 1 : 

sin^ T r 7 e T,\ t>0, (3.10) 

W 7T J r 2 ^ + 2rP cos(/3vr) + 1 ~ ' K J 



o 



and 



sin^vr) 7 r*e ^ t>Q> 







A special case is (3 = \ for which it is known that 

2 1 2 

Ei/ 2 (-Vi) =e t eifc(Vt) = e t ^ / e~ M d«, t>0, (3.12) 
where erfc denotes the complementary error function. 

It may be instructive to note that for sufficiently small times \I/(t) exhibits a 
behaviour similar to that of a stretched exponential; in fact we have 

^(-^) - 1 " f(^T) " ex P{-^/r(l + 13)}, 0<i«l. (3.13) 



Hereafter, we consider two relevant forms for the survival probability ^(t) 
(that we shall denote by fp{t) and gp{t) to distinguish them from ep(t) : = 
Efsi^—t 13 )) which, in exhibiting a decreasing behaviour with a power law decay 
for large times, represent alternative candidates for long-memory processes. 
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The simplest function which meets these requirements is expected to be: 



ftW != i + r(i-fl* ' i£0 " (3 - 14) 



so that 

Tip 



t -> 0, 



sin(7r/3) r(l + /3)' ' (3 15) 

sm^Tr) r(g) 

One can infer from the Eqs (3.7) and (3.15) that e / g(t) and fp{t) practically 
coincide for all t > if (3 is sufficiently small, say for < (3 < 0.25 . For 
greater values of (3 , that are relevant in our subsequent empirical analysis, 
their difference is expected to be appreciable in a wide range of time intervals. 

Another possible choice of major statistical relevance is based on the assump- 
tion that the waiting-time pdf ip(t) may be an extremal, unilateral, stable 
distribution with index (3. In this case the Laplace transforms of ip(t) and 
^(t) read 

^j(s) = exp(-s' 3 ) , fr(s) = 1 - exp (~ g/3 ) ; 0</3<l. (3.16) 

s 

By inversion we obtain for t > 

m = ^0-/3,o (-^) , = 1 - (-^) , 0</3<l, (3.17) 

where 0-/?,o , 0-/3,i denote entire transcendental functions (depending on two 
indices), known as the Wright functions, defined in the complex plane by the 
power series 

oo n 

W=g o „, r(A „ + M) . A>-l, ,ec, zee. (3.18) 

For detailed information on the Wright type functions and their Laplace trans- 
forms the reader may consult e.g. Refs. [17,22,23]. We note that for < (3 < 1 
and /x = 0or/x = l — (3 the corresponding Wright functions appear in the 
fundamental solutions of time-fractional diffusion equations (see e.g. Refs. 
[21,22,24] and references therein). 

Hereafter, like for the Mittag-Leffler-type functions (3.5), we quote for the 
Wright-type functions (3.17) their series expansions and asymptotic represen- 
tations. For the waiting-time pdf we have (see e.g. Ref. [25]), 

( >0, (3,9) 
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and 
where 
A = 



V>(t) ~ Ar a exp (-fir 6 ) , *->o, 

0i/(i-/9) n 1 / 2 



2tt(1-/3) 



a = - 2 ^ , 5 = (1-/3) /3 b , 6 = 



2(1-/5) 



1-/3 



(3.20) 



(3.21) 



For the survival probability we obtain 



71=1 



and 



where 



0/3 



C 



(t) = tf(t) ~ 1 -Ct c exp(-fir 6 ) , i->0, 



2tt(1-/3) /3V(i-/3) 



1/2 







2(1-/3) 2 



(3.22) 

(3.23) 
(3.24) 



Like for the Mittag-Leffler function, a special case is again f3 — \ for which 
we obtain the analytical expressions 



^(t) = -=t 
V 71 " 



-3/2 



exp I - Jt 



#(f) = erf 



2^, 



i > 0. 



(3.25) 



We note that in this particular case the asymptotic representation (3.20) and 
(3.21) provides the sum of the series (3.19) and henceforth the exact expression 
for ip{t) in Eq. (3.25), the so-called Levy-Smirnov pdf (see e.g. Ref. [26]). 

Hereafter we would like to point out the major differences between the Mittag- 
LefHer-type function ep(t) := E^—t 13 ) and the Wright type function gp(t) : = 
l-(f)- f3:1 (-t- f3 ) , that can be inferred by analytical arguments. The first differ- 
ence concerns the decreasing behaviour before the onset of the common power 
law regime: whereas ep(t) starts at t — vertically (the derivative is — oo) and 
is completely monotone, gp(t) starts at t — horizontally (the derivative is 0) 
and then exhibits a change in the concavity from downwards to upwards. A 
second difference concerns the limit for f3 — > 1 ; whereas ep(t) tends to the ex- 
ponential exp(— t) (no memory), gp(t) tends to the box function H(t)—H(t—1) 
(as directly obtained from the Laplace inversion of Eq. (3.16) for (3 = 1). As 
a consequence, the corresponding waiting-time pdf tends to the Dirac delta 
function S(t — 1) , a peculiar case considered by Weiss [9] in his book in p. 47 
as an example of a non-Markovian process. 
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4 Empirical analysis 



In order to corroborate the theory presented above, we have analyzed the 
waiting-time distribution of BUND futures traded at LIFFE in 1997. BUND 
is the German word for bond. Futures are derivative contracts in which a party 
agrees to sell and the other party to buy a fixed amount of an underlying asset 
at a given price and at a future delivery date. In this case the underlying asset 
is a German Government bond. 

We have considered two different delivery dates: June 1997 and September 
1997. Usually, for a future with a certain maturity, transactions begin some 
months before the delivery date. At the beginning, there are few trades a day, 
but closer to the delivery there may be more than 1 000 transactions a day. 
For each maturity, the total number of transaction is greater than 160 000. 

In Figs. 1 and 2 we plot ^(r) for the June and September delivery dates, 
respectively. The circles refer to market data and represent the probability of 
a waiting time greater than the abscissa r. We have determined about 600 
values of ^(r) for r in the interval between 1 and 50 000 s, neglecting the 
intervals of market closure. The solid line is a two-parameter fit obtained by 
using the Mittag-Leffler-type function 

*(r) = e^( 7 r) = Ep [-^rf] , (4.1) 

where (3 is the index of the Mittag-Leffler function and 7 is a time-scale factor, 
depending on the time unit. For the June delivery date we get an index (5 = 
0.96 and a scale factor 7 = ^ , whereas, for the September delivery date, we 
have j3 = 0.95 and 7 = ^ . The fit in Fig. 1 has a reduced chi square x 2 — 0.26, 
whereas the reduced chi square of the fit in Fig. 2 is x 2 ~ 0.25. The chi-square 
values have been computed considering all the values of ^ . 

In Figs. 1 and 2, the dash-dotted line is the stretched exponential function 
exp{ — ( 7 r) /3 )/r(l-|- / 9)} (see Eq. (3.13)), whereas the dashed line is the power- 
law function ( 7 r) _/3 /r(l — (3) (see the second equation in Eq. (3.7)). The 
Mittag-Leffler function interpolates between these two limiting behaviours: 
the stretched exponential for small time intervals, and the power-law for large 
ones. 



Even if the two fits seem to work well, some words of caution are necessary. The 
Mittag-Leffler-type function ep^r) naturally derives from our assumption on 
the "memory function" in the CTRW model; however, as previously observed, 
it is not the unique possibility compatible with a long-memory process with 
a power-law decay. As a consequence, hereafter, we shall also consider for 
^(t) the two alternative functions discussed in Section 3, namely the rational 
function fp{^r) (see Eqs. (3.14) and (3.15)), and the Wright function gp{^r) 
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(see Eqs. (3.22)-(3.24)). 



In Figs. 3 and 4, by taking the same data as in Figs. 1 and 2 respectively, we 
compare the functions ep^r) (solid line), fpi'jr) (dash-dotted line) and gpi'jr) 
(dashed line). Whereas in the previous figures we have adopted a log- log scale 
to point out the power-law decay by a straight-line, now we find it convenient 
to use a linear scale for the ordinates to point out the behaviour of the functions 
for small values of r . From these figures we can infer that the Mittag-Leffler 
function fits the data of the empirical analysis much better than the other two 
chosen functions, thus corroborating our approach to CTRW based on the 
fractional-time derivative. However, the Mittag-Leffler fit significantly differs 
from the empirical data for small values of r . 



5 Conclusions 

The CTRW is a good phenomenological description of the tick-by-tick dynam- 
ics in a financial market. Indeed, the CTRW can naturally take into account 
the pathological time-evolution of financial markets, which is non-Markovian 
and/or non-local. From this point of view, by a proper choice of a (perhaps 
non- stationary) joint pdf r), one could accurately reproduce the statisti- 
cal properties of market evolution. In this respect, the model can be useful for 
applications where Monte-Carlo simulations of market settings are needed. 

With additional assumptions, the CTRW hypothesis can be tested against 
empirical data, thus providing useful information on the restrictions of the 
premises. In this paper, we have assumed a particular form for the time- 
evolution kernel, leading to a time-fractional Kolmogorov-Feller equation. In 
its turn, this implies that ^(t), the probability of finding a waiting-time inter- 
val greater than r, is a Mittag-Leffler function. There is a satisfactory agree- 
ment between this prediction and the empirical distributions analyzed in Figs. 
1-4, but not for small time intervals. 

Among the various questions for future research on this topic, two are partic- 
ularly relevant in our opinion. The first one concerns the behaviour of other 
assets. Prices of liquid stocks should have a completely different time scale. In- 
deed, for the futures here considered, at the beginning of their lifetime several 
hours passed between two consecutive trades, a feature which is not likely to 
be shared by liquid stocks. The second problem concerns the uniqueness of the 
kernel. The Mittag-Leffler kernel yields the elegant time-fractional Kolmogorov 
equation, but there might be other possibilities for interpolating between the 
small waiting-time and the large waiting-time behaviour of \&(t). 

Finally, there is an implication for microscopic market models, a realm where 
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many physicists have started researching. We believe that a microscopic model 
should, at least phenomenologically, take into account that agents in the mar- 
ket decide to sell and buy an asset at randomly distributed instants. It would 
be a success to derive the "right" waiting-time distribution from first princi- 
ples, whatever these first principles will be. 



Appendix. The Caputo fractional derivative 

For the sake of convenience of the reader here we present an introduction to 
the Caputo fractional derivative starting from its representation in the Laplace 
domain and pointing out its difference with respect to the standard Riemann- 
Liouville fractional derivative. So doing,g we avoid the subtleties lying in the 
inversion of fractional integrals. If f(t) is a (sufficiently well-behaved) function 
with Laplace transform L {/(i); s} — f(s) — / °° e f(t) dt , we have 

^{ *} = ^/(*)-^ _1 /(o + ), o</3<i, (Ai) 

if we define 

d 13 

We can also write 

A modification of Eqs. (A.1)-(A.4) holds any non integer /3 > 1 (see [19]). We 
refer to the fractional derivative defined by Eq. (A. 2) as the Caputo fractional 
derivative, since it was formerly applied by Caputo in the late 1960s for mod- 
elling the dissipation effects in Linear Viscoelasticity (see e.g. Refs. [27,22]). 
The reader should observe that this definition differs from the usual one named 
after Riemann and Liouville, which is given by the first term in the R.H.S. 
of (A. 4) see e.g. [28]. For more details on the Caputo fractional derivative we 
refer to Refs. [19,29,30]. 



v 

I 



df(r) dr 



r(l -13) J dr (t- rf ' 



(A2) 
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Fig. 1 

Survival probability for BUND futures with delivery date: June 1997. 
The Mittag-Lefher function (solid line) is compared with the stretched 
exponential (dash-dotted line) and the power (dashed line) functions. 

(0 = 0.96, 7 = 1/12) 




t {s} 



Fig. 2 

Survival probability for BUND futures with delivery date: September 1997. 
The Mittag-Lefner function (solid line) is compared with the stretched 
exponential (dash-dotted line) and the power (dashed line) functions. 

(0 = 0.95, 7 = 1/12) 
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10 10 10 10 10 10 

X {s} 
Fig. 3 

Survival probability for BUND futures with delivery date: June 1997. 
The Mittag-Lemer function (solid line) is compared with the rational 
(dash-dotted line) and the Wright (dashed line) functions. 
(P = 0.96, 7 = 1/12) 




x {s} 
Fig. 4 

Survival probability for BUND futures with delivery date: September 1997. 
The Mittag-Lemer function (solid line) is compared with the rational 
(dash-dotted line) and the Wright (dashed line) functions. 
(^9 = 0.95, 7 = 1/12) 
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